www.gusucode.com > 循环自相关函数工具箱源码程序 > matlab代做 修改 程序循环自相关函数工具箱/cyclostationary_toolbox/ex5.m
clear; clc; %A=1;B=1; fb=7;fa=100;fs=600;N=1024; t=0:1/fs:1024/fs; x=cos(2*pi*fa*t+cos(2*pi*fb*t)); figure(1) plot(t,x); axis([0 1 -1 1]); xlabel('\fontsize{9}\fontname{标楷体}时间/s');ylabel('\fontsize{9}\fontname{标楷体}幅值/um'); x1=x.*x; y=fft(x1,N); Pyy=sqrt(y.*conj(y))/N; figure(2) ff=fs*(0:N/2)/N; plot(ff(1:N/2+1),Pyy(1:N/2+1)); xlabel('\fontsize{9}\fontname{标楷体}循环频率/Hz');ylabel('\fontsize{9}{R_{x}}^{\alpha}(0) /um^2')